SCATTERING BOUNDARY RIGIDITY IN THE PRESENCE OF 

A MAGNETIC FIELD 



PILAR HERREROS 

Abstract. It has been shown in DPSU that, under some additional assump- 
tions, two simple domains with the same scattering data are equivalent. We 
show that the simplicity of a region can be read from the metric in the bound- 
ary and the scattering data. This lets us extend the results in |DPSUI to 
regions with the same scattering data, where only one is known apriori to be 
simple. We will then use this results to resolve a local version of a question 
by Robert Bryant. That is, we show that a surface of constant curvature can 
not be modified in a small region while keeping all the curves of some fixed 
constant geodesic curvatures closed. 



1. Introduction 

A magnetic field on a Riemannian manifold can be represented by a closed 
2 -form 0, or equivalently by the (1,1) tensor Y : TM -> TM denned by fi(£, v) = 
(Y(£),i/) for all x £ M and f, v £ T X M . The trajectory of a charged particle in 
such a magnetic held is then modeled by the equation 

we will call such curves magnetic geodesies. In contrast to regular (or straight) 
geodesies, magnetic geodesies are not reversible, and can't be rescaled, i.e. the 
trajectory depends on the energy \y\ 2 . 

Magnetic geodesies and the magnetic flow where hrst considered by V.I. Arnold 
[Arl] and D.V. Anosov and Y.G. Sinai [AS] . The existence of closed magnetic 
geodesies, and the magnetic flow in general, has been widely studied since then. 
Some of the approaches to this subject are, the Morse-Novikov theory for varia- 
tional functional (e.g. [No ! lNT [ ITa] ) , A ubry-Mather's theory (e.g. [CMP] ), the 
theory of dynamical systems (e.g. [Grog[ lNi| IPPj ) and using methods from sym- 
plectic geometry (e.g. [Ar2l IGiTI lGi2] ) . 



In the special case of surfaces, where the 2— form has the form = k(x)dA, a 
magnetic geodesic of energy c has geodesic curvature k g — k{x)/^fc. This relates 
magnetic geodesies with the problem of prescribing geodesic curvature, in particular 
with the study of curves of constant geodesic curvature. This relation was used by 
V.I. Arnold in |Ar3j . and later by many others (see e.g (Le[ ISchj ). to study the 
existence of closed curves with prescribed geodesic curvature. 

It is clear that on surfaces of constant curvature the curves of large constant 
geodesic curvature are circles, therefore closed. The study of these curves goes back 
to Darboux, who in 1894 claimed (in a footnote in his book |Daj ) that the converse 

l 
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is true, that is, if all curves of constant (sufficiently large) geodesic curvature are 
closed, then the surface has to be of constant Gauss curvature. 

The proof of this result depends strongly on the fact that curves of geodesic 
curvature are closed for all large curvature, or equivalently low energy. This raises 
the following question, brought to my attention by R. Bryant. 

Question 1. Are surfaces of constant Gauss curvature the only surfaces for which 
all curves of a fixed constant nonzero geodesic curvature are closed? 

For the case where the constant is 0, this question corresponds to existence of 
surfaces all of whose geodesies are closed. The first examples of such surfaces where 
given by Zoll [Zo] who, in 1903, constructed a surface of revolution with this prop- 
erty. 

Using the relation between geodesic curvature and magnetic geodesies we can 
approach this question by studying magnetic geodesies on a surface, in the presence 
of a constant magnetic field. A first step in this direction is to determine if a sur- 
face of constant curvature can be locally changed keeping all the curves of a fixed 
constant geodesic curvature closed. We show in section [5] that the metric can not 
be changed in a small region without loosing this property. In fact, more generally, 
we show that a Riemannian surface with a magnetic flow whose orbits are closed 
can't be changed locally without "opening" some of its orbits. 

The easier way of changing the metric in a small region without opening the 
orbits is to require that all magnetic geodesies that enter the region leave it at the 
same place and in the same direction as before, to join the outside part of the orbit. 
This is the magnetic scattering data of the region; for each point and inward direc- 
tion on the boundary, it associates the exit point and direction of the corresponding 
unit speed magnetic geodesic. 

With this problem in mind, we can ask the following boundary rigidity question 
for magnetic geodesies. In a Riemannian manifold with boundary, in the presence 
of a magnetic field, is the metric determined by the metric on the boundary and 
the magnetic scattering data? 

In general this is not true, even for the geodesic case. For example, a round sphere 
with a small disk removed has the same scattering data as a round M.P 2 with a disk 
of the same size removed. One of the usual conditions to obtain boundary rigidity 
is to assume that the region is simple. In our setting simple means a compact 
region that is magnetically convex, and where the magnetic exponential map has 
no conjugate points (see section 2]). 

For simple domains, scattering rigidity for geodesies is equivalent to distance 
boundary rigidity (see |Cr3j ) and it has been widely studied. It is known to hold 
for simple sub-domains of 1™ }Groj or |Cr2j , an open round hemisphere |Mij , hy- 
perbolic space |BCG| }Cr3] , and some spaces of negative curvature |Ot| ICrlj among 
others. For a discussion on the subject, see |Cr3j . 

Recently N. Dairbekov, P. Paternain, P. Stefanov and G. Uhlmann proved (in 
DPSU ) magnetic boundary rigidity between two simple manifolds in several classes 
of metrics, including simple conformal metrics, simple analytic metrics, and all 
2— dimensional simple metrics. 
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FIGURE 1. A closed orbit going through R twice. 

To be able to apply this results to local perturbations of existing metrics we 
would like to be able to compare a simple domain with any other (not necessarily 
simple) domain with the same boundary behavior. For this we prove the following 
theorem. Thus, proving magnetic rigidity for the simple domains considered in 
|I)l'Sl.| . 

Theorem 1. Magnetic simplicity can be read from the metric on the boundary and 
the scattering data. 

To change the metric in a small region without opening the orbits, it is not 
necessary to preserve the scattering data. It could be the case, in principle, that 
orbits exit the region in a different place, but after some time, came back to the 
region and leave it in the proper place to close up again (see figure [T]). In section [5] 
we look at this case in 2 dimensions, and we show the following theorem. 

Theorem 2. Let M and M be compact surfaces with magnetic fields, all of whose 
magnetic geodesies are closed. Let R C M be a strictly magnetically convex region, 
such that every magnetic geodesic passes through R at most once. 

If the metric and magnetic fields of M and M agree outside R, then they have 
the same scattering data. 

Applying these two theorems for a constant magnetic field on a surface of con- 
stant curvature, we conclude that: 

Corollary 3. Given a surface of constant curvature and a fixed k such that all the 
circles of geodesic curvature k are simple. The surface can not be perturbed in a 
small enough region while keeping all the curves of geodesic curvature k closed. 

Acknowledgments. I would like to thank my advisor Christopher Croke for all 
his support. 
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2. Magnetic Jacobi fields and conjugate points 

Let (M, g) be a compact Riemannian manifold with boundary, with a magnetic 
field given by the closed 2-form fi on M. Denote by luo the canonical symplectic 
form on TM, that is, the pull back of the canonical symplectic form of T*M by the 
Riemmanian metric. The geodesic flow can be described as the Hamiltonian flow 
of H w.r.t. ljq, where H : TM — > K is defined as 

H(v) = \\v\% v e TM. 

In a similar way, the magnetic flow ip* : TM — > TM can be described as the Hamil- 
tonian flow of H with respect to the modified symplectic form uj = us Q + 7r*Sl. This 
flow has orbits t —> (7(t), j'it)), where 7 is a magnetic geodesic, i.e. Vy 7' = ^(7'). 
Note that when f2 = we recover the geodesic flow, whose orbits are geodesies. 

It follows from the above definitions that the magnetic geodesies have constant 
speed. In fact, 

| (Y(t), 7 '(t)) = 2 <F( 7 '), 7') = 20(7', 7') = 0. 

Moreover, the trajectories of the magnetic geodesies depend on the energy level. 
Unlike geodesies, a rescaling of a magnetic geodesic is not longer a magnetic ge- 
odesic. We will restrict our attention to a single energy level, or equivalently to 
unit speed magnetic geodesies. Therefore, from now on, we will only consider the 
magnetic flow ijj* : SM SM. 

The choice of energy level is not a restriction, since we can study other energy 
levels by considering the form (l = AS1, for any A £ 1. 

For x € M we define the magnetic exponential map at x to be the partial map 
exp% :T X M M given by 

acp£(tO=W*(0, ^S X M. 

This map takes a vector t£ e T X M to the point in M that corresponds to following 
the magnetic geodesic with initial direction £, a time t. This function is C°° on 
T X M \ {0} but in general only C 1 at 0. The lack of smoothness at the origin can 
be explained by the fact that magnetic geodesies are not reversible. When we pass 
through the origin we change from 7^ to 7_^ , that in general only agree up to first 
order. For a proof see Appendix A in DPSU . 

We will say that a point p S M is conjugate to x along a magnetic geodesic 7 if 
P — 7(^0) — ex Px(to£.) and v = to£ is a critical point of exp%. The multiplicity of 
the conjugate point p is then the dimension of the kernel of d v exp^. 

In what follows, and throughout this paper, if V is a vector field along a geodesic 
7(i), V will denote the covariant derivative V 7 'U. 

We want to give an alternative characterizations of conjugate points. For this 
consider a variation of 7 through magnetic geodesies. That is 



/(*,«) =7.(*) 
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where 7 S (£) is a magnetic geodesic for each s <E (— e, e) and £ £ [0,T]. Therefore 
J| g£ = ^( )■ Using this and the definition of the curvature tensor we can write: 

^fW^ - D_D_dl _ D_D_d£ _ R( d£ df\d£ 

d s y 1 yat)) — ds at at at as at „ ^ds* at' at 

— S.S.2L 4- U(2L SL\<1L 

— at at as ^ n \ at > as J at 

If we call the variational field J(t) = §f (i, 0) we get for s = 

v J (y( 7 / )) = ^" + ^(7 , J «/)7'- 

Note also that 

Vj(Y( 7 ')) = Y(VjY) + (VjYW) 

and 

so we can rewrite the above equation as 

J" + R(j', J)i - y(j') - (v J y)( 7 / ) = o. 

Since magnetic geodesies can't be rescaled we have |7g| = 1 for all the magnetic 



geodesies in the variation. This equation together with equation (|2.ip gives, for any 
such variational field J, 

(J',Y) = (Vj'y', 1 ')=J(Y,Y) = 0. 

This equations characterize the variational field of variations through magnetic 
geodesies, in a way analogous to the characterization of Jacobi fields. We will use 
this equation as a definition as follows. 

Given a magnetic geodesic 7 : [0, T] —> M, let A and C be the operators on 
smooth vector fields Z along 7 defined by 

A(Z) = Z" + R(j', Z)i - Y(Z') - (V Z Y)( 7 % 

C{Z) = R{ 1 ',Z)i-Y{Z')-{V z Y){i). 

A vector field J along 7 is said to be a magnetic Jacobi field if it satisfies the 
equations 

(2.2) A(J) = 
and 

(2.3) (J',7')=0. 

Note that from equation 12.21 we can see that 

&{JW) ={J",i) + {J',Y{i)) 

= (-i?( 7 ', J)i + Y(J') + (VjF)(Y), 7') - (Y(J'), 7') = 

where we used that ((V jY)(-f'), 7') = because Y is skew-symmetric. Therefore, 
it is enough to check condition 12.31 at a point. 

Note: In [DPSU] magnetic Jacobi fields were defined without condition 12.31 al- 
though where used only for vector fields that satisfied this condition. We will add 
it to the definition since it helps preserve the relation between (magnetic) Jacobi 
fields and variations through (magnetic) geodesies. 
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A magnetic Jacobi field along a magnetic geodesic 7 is uniquely determined by 
its initial conditions J(0) and J'(0)- To see this, consider the orthonormal basis 
defined by extending an orthonormal basis e±, . . . , e„ at 7(0) by requiring that 

(2.4) ei = Y(ei) 

along 7. This extension gives an orthonormal basis at each point since 
±{e i ,e j ) = {Y(e i ),e j ) + {e i ,Y(e j ))=0. 

Using this basis, 

n 

Z = ^2 fi e i 
i=l 

and we can write equation 12.21 as the system 

n n 
i=l i=l 

where yij = (Y(ei),ej) and 

ay = (VyYfe) + i?( 7 ', eijV - Y(Y(e t )) - (V^YXY), e,-) . 

This is a linear second order system, and therefore it has a unique solution for each 
set of initial conditions. 

Magnetic Jacobi fields correspond exactly to variational field of variations through 
magnetic geodesies. In the case of magnetic Jacobi fields J along 7^ that vanish at 
this can be seen by considering 

(2.5) /(*,*) =7.(*) = 

where f : (-e,e) ->■ S X M is a curve with £(Q) = f = 7c and £'(Q) = J'(0). This 
is clearly a variation through magnetic geodesies, and therefore its variational field 
§f(i, 0) satisfies 12.21 The variational field 0) and the magnetic Jacobi field 
J(t) are then solutions of 12.21 with the same initial conditions, therefore they must 
agree. 

For magnetic Jacobi fields J that do not vanish at 0, we can use the variation 

/(i, a) = 7 .(i) = ezp£ w (<£(«)) 

where t(s) is any curve with r'(0) = J(0) and £(s) is a vector field along r with 
C(O) = Y(0) and ^'(0) = J'(0). 

It is easy to see from the definition and the equation 7" = 1^(7') that 7' is 
always a magnetic Jacobi field. Unlike the case of straight geodesies, this is the 
only magnetic Jacobi field parallel to 7'. Another difference from the straight 
geodesic case is that magnetic Jacobi fields that are perpendicular to 7' at t = 
don't stay perpendicular for all t. For this reason we will sometimes consider instead 
the orthogonal projection J- 1 = J — fj' where / = (J, 7'). The component fj' of 
J parallel to 7' is uniquely determined by and J(O), since 

/' = (J', 7') + (J, 7") = (J,Y( 7 ')) = (J X ,Y( 7 ')) . 
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We will need one more property of magnetic Jacobi fields that vanish at 0. Let 
7 : [0, T] — >• M be a magnetic geodesic with 7(0) = x. Let v £ T y /S X M, or 
equivalently under the usual identification, v £ T X M perpendicular to 7'. 

Let f(t, s) be a variation through magnetic geodesies of the form 12.51 where 
t G [0,T] and £ : (-e,c) -> S X M with £(0) = Y(0) and f (0) = v. The variational 
field J.„ of this variation is 

J vH) =f (*.0)= &[W t K(*))]|^ = WW t ](£ / (0)) 

(2.6) 

= d^t^n o d^i)r{v) = <Ly>(t)ir ° rf 7 '(o)^ («) 
and its derivative is given by 

■£(*) =I|f(*,0)=£|[ 7 ro^(e( S ))]| s=0 

= £[^(^( s ))]| s =o = ^(0)^(^(0)) = d y( o)^(«). 
This equations are independent of the variation /. 

Since the magnetic flow is a Hamiltonian flow with respect to the symplectic 
form u> = luq + 7r*f2, this form is invariant under the magnetic flow ijj* [Pal pg. 10]. 
Therefore for any two magnetic Jacobi fields J v and J w as above, we have that 

u}(<L y r^ t (v),<Ly>ip i (ui)) = U3o(<L,r^ t (v), d^'ip^w)) + %*Q,((L 1 >ip t {v), d 1 >^ t {w)) 

= (j v (t), j' w (t)) - (.r v {t), j w {t)) + n(j v (t), j w (t)) 

is independent of t. Using also that J v (0) — 0, we get 
(2.7) (J„, J' w ) - {J' v , J w ) + (Y(J V ), J w ) = 

for any two such Jacobi fields. 

We will now relate the concepts of magnetic Jacobi fields and conjugate points. 

Proposition 2.1. Let 7^ : [0, T] — > M be the magnetic geodesic with 7(0) = x and 
7'(0) = £. The point p = 7 (to) * s conjugate to x along 7 if and only if there exist 
a magnetic Jacobi field J along 7, not identically zero, with J(0) = and J(io) 
parallel to 7'. 

Moreover, the multiplicity of p as a conjugate point is equal to the number of 
linearly independent such Jacobi fields. 

Consider a variation through magnetic geodesies as in (|2.5[) . with £'(0) = v 
perpendicular to 7'. Then 

df 

Mt) = -±(t,Q) = dttexp£(tv) 

OS 

is a nontrivial magnetic Jacobi field. If there is a vector v for which J v (to) is parallel 
to 7', then dt £exp%(tov) and dt^expf^) = 7' will be parallel, and to£ is a critical 
point of expif.. Conversely, if to£ is a critical point there must be a vector v such 
that d to ^expf.{v) — 0. Let v 1 - = v — (v, £) £, this is not since dt^expf.^) — 7' 7^ 0, 
and 

J„x(i ) = d toi exp^(t v) - d toi exp^(t £) = -t (t>,£) 7' (to). 

To prove the second statement, note that Jacobi fields J Vi as above are linearly 
independent iff the vectors vt are. Since all Vi are perpendicular to 7'(0), the num- 
ber of linearly independent vectors will be the dimension of the kernel of dt ^exp^, 
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that is the multiplicity of the conjugate point. 



3. The Index form 

Let A denote the K- vector space of piecewise smooth vector fields Z along 7. 
Define the quadratic form Ind : A — > R by 

Ind^Z) = ( T {\Z'\ 2 - (C(Z),Z) - (Y(j'),Z) 2 }dt. 
Jo 

Note that 

Ind^Z) = - f T {(A(Z),Z) + (Ytf), Zf}dt + (Z, Z% + V (z, z'- -z'+)\ t ,. 
Jo 

where Z r± stands for the left and right derivatives of Z at the points ti where the 
derivative is discontinuous. 

The Ind 1 {Z) generalizes the index form of a geodesic in a Riemannian manifold. 
It is easy to see that when O = these are the same form. We will see throughout 
this section that, when restricted to orthogonal vector fields, we retain some of the 
relations between (magnetic) Jacobi fields, index form and conjugate points. 

Let Ao denote the R-vector space of piecewise smooth vector fields Z along 7 
such that Z(0) — Z(T) = 0, A 1 - the subspace of piecewise smooth vector fields that 
stay orthogonal to 7', and Aq = Aq n A^. 

For any magnetic Jacobi field J along a magnetic geodesic 7, let / = (J, 7') and 
J 1 - = J — / 7 ' the component of J orthogonal to 7 , using that 7" = F( 7 ') we have 
A(J) = 

J ± %/V+2/'r(70+/7'''+i?(7',^+/7')7'-^(J ± '+/Y+/7'')-(Vz + / r r)(7') 
so 

(3.1) = A(J^) + fA{i) + f'i + f'Y{i) 

On the other hand, any magnetic Jacobi field satisfies (J', 7') = 0. Using this 
together with / = (J, 7') and 7" = F( 7 ') we see that /' = (J, Y{i)) = ( J x , F( 7 ')>. 
Since 7' is a Jacobi field, A(j') = 0, and we have from (13. ip 

(A(j^z) = -f (y(V), j±) = (r( 7 0, J^) 2 

therefore if J is a magnetic Jacobi field and its orthogonal component J 1 - is in Ao 

(3.2) Ind 1 (Z) = - {(A(J ± ),J ± ) + (Y( 1 '),J ± ) 2 }dt = 0. 

Jo 

Also, if there are no conjugate points along 7, we can restrict our attention to 
a neighborhood of 7 for which the magnetic exponential is a diffeomorphism. In 
this case we can adapt the proof given in [DPSU for the case of simple domains to 
prove the following version of the Index Lemma. We include the proof below. 
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Lemma 3.1. Let 7 be a magnetic geodesic without conjugate points to 7(0). Let J 
be a magnetic Jacobi field, and J 1 - the component orthogonal to 7'. Let Z £ A 1 - be 
a piecewise differentiate vector field along 7, perpendicular to 7'. Suppose that 

(3.3) J^O) = Z(0) = and J^{T) = Z(T) 

Then 

Ind~,{J^) < Indy(Z) 
and equality occurs if and only if Z = J 1 - . 

Note that in the case where the vector field Z satisfies Z(0) = Z(T) = 0, J = 7' 
is a Jacobi field that satisfies the above hypothesis, giving the following corollary. 

Corollary 3.2. If 7 has no conjugate points and Z £ Aq, then Ind 1 (Z) > 0, with 
equality if and only if Z = 0. 

In other words, if 7 has no conjugate points, the quadratic form Ind 1 : Aq — ► K 
is positive definite. 

Proof of Lemma \3.1\ Given a vector v £ T 7 ( )M orthogonal to 7', we can define 
a magnetic Jacobi field J v along 7 as in (|2.6|) that has J v (0) = and J^,(0) = u. 
Since there are no conjugate points to 7(0) along 7, these Jacobi fields are never 
nor parallel to 7'. There exist a basis {v\, . . . ,v n } of T^tmM, with v\ — 7' such 
that J\(t) — "/(£) and Jj(i) = J Vi (t), i = 2,. . . ,n form a basis for T^n\M for all 

ie(o,T]. 

If Z is a vector field with Z(0) — 0, we can write, for t G (0,T] 

n 

z(*) = X)/i(*) J i(*) 

i=l 

where /i,...,/ n are smooth functions. This function can be smoothly extended 
to t = 0, as we now show. For i > 2, we can write Jj(t) = £vlj(t) where A4 are 
smooth vector fields with Ai(0) — Jj'(0). Then each is parallel to Ji(t) and 

{y(t), A 2 (t), • • .,A n (t)} is a basis for all t £ [0,T], so 

n 

Z(t) = gil ' + J29i(t)Ai(t), 

It follows that for t £ (0, T], gi = fx and g t (t) = tfi(t) for i > 2, and since Z(0) = 0, 
<7t(0) = and /, extends smoothly to t — 0. 

Using this representation we can write 

Ind~f(Z) = 

e r (- 4 (^^)' /^i) r ^> 2 *+ (z, z% ( z > z '~ z ' + )\ th ■ 

where Z ,=t stands for the left and right derivatives of Z at the points tk where the 
derivative is discontinuous. 
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And 
(3.4) 

AtfiJi) = fl'Ji + 2f!j! + UJ'J + fcRW, Ji)i ~ flY(Ji) - fiY(Ji) - h{V . h Y){i) 
= fiA(J i ) + fl'J i + 2ftJ>-f!Y(J i ) 

where the first term is since Jj is a magnetic Jacobi field. 



We know, moreover, from (|2 . 7[) that 

(J i ,J' j )-(Ji,J j ) + (Y(J i ),J j )=0 

so we can write 

(A(fiJi), Jj) = (fi'Ji, J,) + 2 (f<Jl, J,) - (flY(Ji), Jj) 
= ti' (Ji,Jj) + fl(Jl,Jj) + fl(Ji,JS) 

and 

T (A(fiJi), fjjj) dt = £ f j t (f t (J h J 3 ))dt = (flJi, fjJj^-f (ftJi, f'jJj) dt. 
Using this, and that Z(0) = we can write the index form as 

T In 







Ind. 



( z )= / nE/^ii 2rf *- (E^> z 

JO ! \ ! 

- [ T (Y(Y),zf dt + (Z, Z')\ T Q + ( Z > Z " z ' + ) 
Jo 

„T n In \ I n 1 

= j o iiE/ l '^ii 2 -( r (70^> 2 *+^(r)E/ 4 ( r )^KT)UE(^E^( j r-^ 

where the last term is because Jj is diffcrcntiable. 

Let W — f'iJii an d remember that J\ = 7', then we get 

n 

II E f'^W 2 = <#V + w > /iV + WO = ^ 2 + 2 ^ (V. wo + (W, W) . 

1 

Since Z is orthogonal to 7', (Z,j') = and by differentiating (Z',-f') — — (Z,Y(j')). 
Using also that for magnetic Jacobi fields (J', 7') = we have 

- (Z, F( 7 ')) = ^E ^ + A = fi + (W, 7) 

and 



so 

{•T n 



(Z,Y{ 1 ')) l =f[ 2 + 2f[(W r i') + (W n ') 2 . 
r-r n / n \ 

/nd 7 (z) = ^ ||^//j i || 2 -<y(y),z) 3 * + /z(nE^ T ) J *( T )) 

j\w,W)-(W^fdt+/z(T),j2n(T)J^T)\ 
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(3.5) 



£ ' 11^11*+ /z(T),j2Mnn(T) 



where W 1 - is the component of W orthogonal to 7'. 

Let J be the magnetic Jacobi field J — h{T)Ji- The fact that it is a magnetic 
Jacobi field is easy to see since A is linear over M., it is also clear that J(0) = and 
J-'-(T) = Z{T). Since 

n n 

2 2 

the corresponding functions fi(T) are constant for i > 2, so W — and we can see 
from equation (|3.5p that 

Ind y (J ± ) = (j(T),^f i (T)Ji{T)\ 

that gives 



Ind 7 (Z) - Ind^J^) = \\W ± \\dt>0 
Jo 

with equality iff vanishes everywhere. That is 

n 

2 

and therefore fi constant for i > 2. So Z = /17' + /j J,; = /17' + J, and since 
Z is orthogonal to 7' this implies that Z = J . 

□ 

In what follows we will want to use the above lemma for more general vector 
fields, that do not vanish at but vanish at T. Since magnetic flows are not 
reversible, we can't simply reverse time. We will consider instead the associated 
magnetic flow (M, g, — ft). This magnetic flow has the same magnetic geodesies, 
but with opposite orientation. 

Lemma 3.3. Let (M, g, Q) be a magnetic field. Then (M, g, — fl) is also a magnetic 
field. 

(1) Magnetic geodesies in both magnetic fields agree, but with opposite orienta- 
tion. 

(2) Jacobi fields agree in both magnetic fields. 

(3) Index form is independent of its orientation. 

If 7 : [0, T] — > M is a magnetic geodesic in (M, g, fi), denote by 7_ the geodesic 
with opposite orientation, that is 7_(t) = j(—t), for t € [— T, 0]. Then 7"(t) = 
7"(— t) = Y"(7'(— t)) = — Y(7^_(i)), so 7_ is a magnetic geodesic in (M,g,— fi). 
Part 2 follows from 1 and the fact that magnetic Jacobi fields are variational fields 
of variations through magnetic geodesies. Alternatively, we can check that for 
J_(i) = J(-t): 

A-(J-(t)) = J'l(t) + Rb-(t),J-(t)h'-(t)+Y{JL(t)) + (Vj_Y)(7i(i)) 

= J"(-t) + fl(Y(-t), J(-t))V(-t) - Y{J\-t)) (VjY)(V(-t)) 

=^(J(-t)). 
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(j:(*),y_(*)> = (-/(-*), -y(-t)) = (j'(-t),7'(-i)> 



From the above computation follows also that C_(Z_(£)) = C(Z(t)). So 
Ind,_(Z-) = /° {|^(t)| 2 -(C_(Z_(t)),Z_(t)>-(y_( 7 UO),Z_(t)) 2 }^ 

J — T 

° {\Z'(-t)f - (C(Z(-t)), Z(-t)> - (F( 7 '(-t)), Z(-t)> 2 }dt 

-T 

T {|Z'(t)| 2 - (C(Z(t)), Z(t)> - (y(y(t)), Z(t)) 2 }dt = Ind 7 (Z). 





Corollary 3.4. Lemma \3.1\ holds when we replace equation (|3.3[) im£/i Z(0) — 
J-L(O) and Z(T) = J X (T) = 0. 

When Z(T) = J ± (T) — 0, we can consider Z_ and J^, this will satisfy the 
hypothesis of lemma [3"7T1 so we have: 

/ratty (J- 1 ) = Ind y _(J±) < Ind 7 _(Z-) = Ind 7 (Z). 



Lemma 3.5. Ifj(to) is conjugate to 7(0) along 7, for some to < T, then there is 
a vector field Z £ Ag with Ind 1 (Z) < 0. 

Let J a Jacobi field along 7 with J(0) = J (to) = 0, and J be J- 1 for i e [0, to] and 

for t £ [to,T]. Then Ind 1 (J) — 0. We can use the Index Lemma and corollary 

13.41 to show that cutting the corner at to by replacing J by a Jacobi held 

[t -e,t +e] 

with the same endpoints decreases the value of the Index form. So this new vector 
held has Ind n (Z) < 0. 

Lemma 3.6. If 7 (T) is the first conjugate point to 7(0) along 7 and Z £ Aq , 
then Ind 1 (Z) > 0, with equality if and only if Z = or Z is the perpendicular 
component of a Jacobi field. 

This is an extension of lemma [3~T1 and it is proved by a similar argument. We 
will follow the proof of lemma [3TT1 using the same notation. Suppose 7 (T) is a first 
conjugate point and has multiplicity k, then we can find a basis {vq, ■ ■ ■ ,w„_i} of 
T 7 ( )M, with w = 7' such that J Vi (T) are parallel to 7 '(T) for i = 1, . . . , k and are 
not parallel to j'(T) for i = k + 1, . . . , n- 1. Then J (t) = j'(t) and Ji(i) = J Vi (t), 

1 = 1, . . . , n — 1 form a basis for T^u-\M for all £ € (0, T). 

If Z is a vector field in Aq- , we can write, for t £ (0, T) 

n-l 

Z(t) = 

i=0 

where fo, ■ ■ ■ , /«,— 1 are smooth functions. We can extend these functions to t = as 
before. To extend fo to t = T we can write Jj = (t— T)Ai + Ji(T), for i = 1, . . . , fc. 
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Then Ai are smooth vector fields with Ai(T) = J[{T) that is orthogonal to 7' since 
(J', 7') = for all Jacobi fields. Then {7', A\, . . . , Ak, Jk+i, • ■ • , Jn-i} are a basis 
for all t G [0,T], and 

k n — 1 

i=l i=k+l 



It follows that for t £ (0,T), g { = ft for i > k, g t {t) = (t- T)ft(t) for < i < k, 
and 50 = f + £* =1 (JiCO.Y)- Since Z(T) = 0, ^(T) = and ft extends 
smoothly to t = T. 

Following the proof of lemma 13.11 we get from (|3.5[) 

IndyiZ) = jT \\W x \\dt + (z(T),j^ft(T)Jl(T)\ 

r T 

= / ||^||dt > 

with equality iff W x vanishes everywhere. That is when ft constant for i > 0. So 
Z = /17' + J for some Jacobi field J, and since Z is orthogonal to 7' this implies 
that Z = J ± . 



Corollary 3.7. Ind 1 {Z) restricted to is positive definite if and only if 7 has 
no conjugate points. 

When 7 has no conjugate points, it follows directly from corollary 13 . 2 1 that Ind 1 
is positive definite. In the case that the endpoints are conjugate to each other lnd 1 
has nontrivial kernel, as can be seen from equation 13.21 If 7 has conjugate points, 
we saw on lemma [3~5l that there is a vector field in K§ with Ind-y < 0, therefore it 
is not positive definite. 

We will be interested in the dependence of the index form on its parameters. 
For this consider a continuous (possibly constant) family of vectors £(s) £ S X M 
and the correspondent family of magnetic geodesies J s (t) — exp%(t£(s)). Let T s , 
the length of each geodesic, be continuous on s. Let A s denote the vector space 
of piecewise smooth vector fields Z s along 7 S , perpendicular to 7^ and such that 
Z(0) = Z{T S ) = 0. 

Let {vi, . . . ,v n } be an orthonormal basis with v\ = 7o(0), and extend it to a 
continuous family {vi(s), . . . , v n (s)} of orthonormal basis for each s with vi(s) = 
£(s). This can be done by defining 

Vi(s) = p s {Vi) 

where p s is a rotation of S n with p s (£(0)) = £(s). We extend this for all t by 
requiring that 

(3.6) V 7 ,e 4 = r( ei ) 

along each magnetic geodesic. As in 12.41 this gives an orthonormal basis for each 
point. 
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Using this basis, we can extend any vector field Z = a i(t) e i(t) m Ao to a 
vector field over the family of geodesies by 

2 s 

that belongs to A s when restricted to each 7 S . We will denote the set of such vector 
fields by A [0 ,i]. 

Since everything depends continuously on s, so does 

Ind Js (Z) = f TS {\Zf-(C(Z),Z) - {Y(i s ),Z) 2 }dt. 
Jo 

We will be mostly interested on whether the Index form is positive definite. 
For such a family of curves, the fact that the Index form is positive definite (and 
therefore the non existence of conjugate points) depends continuously on s in the 
following sense. If the index form is positive definite for some so, and has a negative 
value for some s\ there must be some s € (so,s±) where it has non-trivial kernel. 
Moreover, the first such s will be when 7,,. has conjugate endpoints and no conjugate 
points in the interior. 



4. Simple metrics and boundary data 

Consider a manifold Mi such that M C int(Mi), extend g and f2 smoothly. We 
say that M is magnetic convex at x £ dM if there is a neighborhood U of x in Mi 
such that all unit speed magnetic geodesies in U, passing through x and tangent to 
dM at x, lie in U\int(M). It is not hard to see that this definition depend neither 
on the choice of M x nor on the way we extend g and Vt, to Mi. 

Let II stand for the second fundamental form of dM and v(x) for the inward 
pointing normal. Then if M is magnetic convex 

for all (x,£) e TM. [DPSUI Lemma A.6]. 

We say that dM is strictly magnetic convex if 

n(x,o > (y*(0,K*)) 

for all (x,£) £ TM. 

This condition implies that the tangent geodesies do not intersect M except for 
x, as shown in [DPSUi Lemma A.6]. 

We say that M is simple(w.r.t.(g, f2)) if dM is strictly magnetic convex and the 
magnetic exponential map exp% : {exp^.)~ l (M) — > M is a diffcomorphism for every 
x e M. 

For £ SM, let 7^ : [l~(x, £), l(x, £)] — > M be the magnetic geodesic such 

that 7 5 (0) = x, 7^(0) = £, and ^{l {x, £)), ^(l(x, £)) G <9M. Where /- and Z can 
take the values ±00 if the magnetic geodesic 75 stays in the interior of M for all 
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time in the corresponding direction. 



Let d+SM and d-SM denote the bundles of inward and outward unit vectors 
over dM: 

d+SM = {{x,0 E SM :xE dM, (£,v(x)) > 0}, 
dSM = {(>,£) e SM : x e dM, (£,v(x)) < 0}, 

where v is the inward unit normal to dM. Note that d(SM) = d + SM U d-SM 

and d+SM n d- SM = S(dM) . 

In the case that M is simple, is clear that the functions l~(x,£) and l(x, £) are 
continuous and, on using the implicit function theorem, they are easily seen to 
be smooth near a point (x, £) such that the magnetic geodesic 7g(t) meets dM 
transversely at t = l~(x,£) and t = l(x,£) respectively. By the definition of strict 
magnetic convexity, 7^(t) meets dM transversely for all (x,£) € SM\ S(dM). In 
fact, these functions are smooth everywhere, as was shown by Dairbekov, Pater- 
nain, Stefanov and Uhlmann in the following lemma. 



Lemma 4.1. |DPSU1 Lemma 2.3] For a simple magnetic system, the function 
L : d(SM) R, defined by 



h(x,0 := 



l(x,{) if(x,0ed+SM 
J-(x,0 if(x,t)€d-SM 

is smooth. In Particular, I : d+SM —> R is smooth. The ratio 
is uniformly bounded on d(SM) \ S(dM). 



This lemma was proved as stated, for simple magnetic systems, but the proof is 
a local argument using only the strong magnetic convexity of the region. 

The scattering relation S : d+SM — > d-SM of a magnetic system (M,g,Q) is 
defined as follows: 

s(x,o = hdK^O),i'd l (^0)) 

when the value Z(x, £) is finite, otherwise it is not defined. 

The restricted scattering relation s : d+SM — > dM is defined to be the post- 
composition of the scattering relation with the natural projection of d-SM to dM, 
i.e., 

when properly defined. 

We are interested only in simple domains, and domains that have the same scat- 
tering data as a simple domain, so we will assume that l~ and I are finite and 
smooth on d{SM). Moreover, it follows from the smoothness of I and their defini- 
tions that both S and s are smooth everywhere on d+SM . 
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Figure 2. The Scattering relation. 

Let M be a compact simple domain with respect to ft in the interior of a manifold 
(Mi,g). Let M be a compact domain in the interior of a manifold (Mi, g) with ft. 
Related in such way that g — g and ft — ft on M\\M — Mi\M, and the (restricted) 
scattering relations S, S agree on dM(= dM). To be able to compare the magnetic 
flows in M and M we would like to say that M is also simple, without having 
to impose it as a condition. The purpose of this section is to prove the following 
theorem. 

Theorem 4.2. Given M and M as above, then M is also simple. 

To prove that the magnetic exponential is a diffeomorphism we need to show 
that it has no conjugate points. For this we need the following lemma. 



Lemma 4.3. If there are conjugate points in M, then there is a pair of points in 
dM conjugate to each other. 

Suppose there is a point in the interior conjugate to a; € dM along a geodesic 
7g. Let t : [0, 1] — » S X M be a curve joining S x dM to £, and consider the family 
of magnetic geodesies 7 S = exp£(tr(s)). These geodesies exit M at time l(x,Y s ), 
that by the simplicity of M is a continuous function of s. Close enough to x the 
magnetic exponential is a diffeomorphism, and by lemma 14.11 there is a C > such 
that l{x,rj) < C(i/(x),rj) for all n 6 S X M. This implies that for s small enough, 
the magnetic geodesic from x to s(r(s)) is short, and stays inside a neighborhood 
where the magnetic exponential is a diffeomorphism. Therefore, it has no conjugate 
points, and the index form is positive definite close to x. On the other hand, there 
is a perpendicular vector field along 7^ for which the index form is negative. Then 
Ind la is positive definite for s = and not for s = 1. Let sq be the smallest s for 
which Ind la has non trivial kernel. Then, by the results on the previous section, 
s(t(so)) is conjugate to x along the magnetic geodesic 7 S that joins them. 

If there are points conjugate to each other along a magnetic geodesic 7g, and 
both lie in the interior of M, there must be a point conjugate to 7^(0) along this 
magnetic geodesic. Therefore reducing the problem to the case above. This can be 
proved by a similar argument using the family of geodesies 7^ | r Q sT j . 

Proof of Theorem \4-£\ It is easy to see from the definition that the domain M has 
to be strictly magnetic convex, since the metrics and magnetic flows agree outside 
M. 
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To prove that the magnetic exponential map is a diffeomorphism form (exp^)~ 1 (M) 
to M we need to show that it has no conjugate points, i.e. there are no points in M 
that are conjugate to each other along a magnetic geodesic. For this purpose assume 
such points exist, then by Lemma 14.31 there are points dM conjugate to each 

other along a magnetic geodesic 7g, where 7^(0) = x and 7f(io) = V for some to > 0. 

Let J be a magnetic Jacobi field along 7^ that vanishes at 0, and f(s,t) a 
variation through magnetic geodesies with /(0, t) — 7^ and J as a variational field. 
We can use /(s, t) = j s (t) = exp%(t£(s)) where £ : (— e, e) — > S X M is a curve with 
£(0) = f, £'(0) = J'(0). / is well defined in M for (s,t) £ (-e,e) x [0,T S ] where 
T s = /(x,£(s)). Consider c(s) = f(s,T s ) £ 9M the curve of the exit points in dM. 
Then 

dc , , df , , c?f dT, , 

= J(T s ) + ^(0,T s )^. 

If 7£(2(x, £)) is conjugate to 7^(0) along 7^, there is a Jacobi field J that is at 
t = and parallel to 7^ at T s , then (0) is parallel to 7^. On the other hand, if 

(0) is parallel to 7^ for any Jacobi field with J(0) = then J(T S ) is parallel to 
7g. Therefore 7^(T S ) is conjugate to 75 (0) along 7^. 

Note that we can write c(s) = s(£(s)), that depends only on the scattering data, 
so the scattering relation detects conjugate points in the boundary. Since there are 
no conjugate points in the boundary of M, there can be none in M. Therefore the 
magnetic exponential is a local diffeomorphism. 

We will now see that exp% is a global diffeomorphism from (exp^)~ 1 (M) to M. 
To see that it is surjective let x £ dM, and y any point in M . Let c : [0, 1] — > M 
be a path from x to y, and consider the set A C [0, 1] of points such that c(s) is 
in the image of exp%. This set is open, since exp^ is a local diffeomorphism. To 
see that it is closed, choose a sequence s„ e A converging to sq. Then c(s n ) = 
exp^(t(s„)£(s„)), and there is a subsequence such that t(s n ) and £(s„) converge to 
t and £ respectively. If to£o ^ ( ea; ? 5 x)~ 1 (-^)! there must be a first t\ < t such 
that exp^(ti£o) G 5M. Then exp%(t£ ) must be tangent to 9M and inside M for 
t < ti, which contradicts the magnetic convexity of M. Then, A is both open and 
closed, therefore A = [0, 1] and y is in the image of exp^,. 

To see that exp% is injective for x £ 9M, note that it is a covering map. The 
point x has only one preimage, since by the simplicity of M there are no magnetic 
geodesies form x to x. Therefore exp% is a covering map of degree 1. 

To prove this for x ^ dM, we need to see that there are no trapped magnetic 
geodesies, that is, that there are no magnetic geodesies that stay inside M for an 
infinite time. Note that, since any magnetic geodesic that enters the region at 
£ has to exit at s(£), it is enough to see that all geodesies enter the region at a 
finite time. Let 7 be a magnetic geodesic. We know that we can reach the point 
7(0) from the boundary, so there is a variation through magnetic geodesies 7 S (£) 
with 70 = 7, 7 S (0) = 7(0) for all s £ [0,1], and 71 (ti) £ dM for some t x < 0. 
If 7 So intersects dM, by the magnetic convexity of M it has to be a transverse 
intersection, therefore intersecting dM is an open condition on [0,1]. It is also 
a closed condition, by continuity of the geodesic flow and compactness of dM. 
Therefore, since 71 intersects dM, so does 7^ for all s, and 7 is not trapped. 
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Now we see that expf. is a global diffeomorphism from (exp£)~ 1 (M) to M for 
x £ dM. Since expf. is injective for x € 9M, and all geodesies come from some 
point x in dM , magnetic geodesies in M have no self intersections. In particular, any 
x G M has only one preimage under exp'f.. We can then follow the same argument 
as for x G <9M to show that expf is a global diffeomorphism from (exp^) _1 (M) to 
M, for all ieAf. □ 

5. Rigidity for Surfaces 

Consider a magnetic field on a surface M all of whose orbits are closed, and 
consider a magnetically simple region R on it. We want to prove that there is no 
way of changing the metric and magnetic field in this region in such a way that all 
orbits are still closed. 

In the previous section we saw that such a region is magnetically rigid, therefore 
it can't be changed on the region preserving the scattering data. Here we will look 
at the general behavior of such a magnetic flow to ensure that there are no other 
metrics with all its orbits closed. We want to rule out the case where a magnetic 
geodesic that passes through the region, after coming out at a different spot and 
following the corresponding orbit, goes back into the region and exits at the exit 
point and direction of the original first magnetic geodesic, therefore forming a closed 
orbit out of two (or more) segments of the original orbits, like in figure [TJ 

To show this, assume that we have such a magnetic field. Assume, moreover, 
that the region R is such that every magnetic geodesic passes through R at most 
once. This condition restricts both the size of R and the flow, since in a flow where 
the orbits have many self-intersections such a region might not exist. On the other 
hand, if all orbits are simple we can see by compactness that there are small regions 
with this property. 

Proof of Theorem^ Consider the unit tangent bundle SM, and the magnetic geo- 
desic vector field G i.e. the vector field that generates the magnetic flow on the unit 
tangent bundle. For the sake of simplicity of the exposition we will assume first 
that SM is oriented. Then SM is a compact orientable 3-dimensional manifold, 
and G is a smooth vector field that foliates SM by circles. By a theorem of Epstein 
|Ep| , this foliation is C°° diffeomorphic to a Seifert fibration. In particular, any 
orbit has a neighborhood diffeomorphic to a standard fibered torus. 

Note that to each orbit on SM we can uniquely associate a magnetic geodesic, 
by projecting the orbit back to M. We will use this correspondence freely. As a 
Seifert fibration, the base B or space of orbits of SM is a 2-dimensional orbifold. 

Let SR be the subset of SM that corresponds to the region R, and SdR the 
subset of SM corresponding to vectors tangent to the boundary of R. The orbit of 
a point in SdR corresponds to a magnetic geodesic that is tangent to R, and since 
R is strictly magnetically convex, it is tangent only at one point. This magnetic ge- 
odesic corresponds exactly to one in M, and therefore stays away from R thereafter. 
This means that each orbit contains at most one point of SdR, so the set of orbits 
passing through it forms a 1-dimensional submanifold on B, we will denote it by i?o- 

Let m : B — > N be a function that counts the number of times the orbit pases 
through SR in a common period. For regular orbits this is the number of times 
it passes through SR. If the orbit is singular it has a neighborhood diffeomorphic 
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to an (a, b) torus, that is a torus obtained by gluing two faces of a cylinder with a 
rotation by an angle of 2-Kb /a. In this case the common period is a times the period 
of the singular orbit. Therefore, m will be a times the number of times the orbit 
passes through SR. Since the other orbits in the neighborhood will be completed 
when the singular orbit is traveled a times, m will be, in general, continuous at 
such points. In fact, if a magnetic geodesic intersects dR transversally (or not at 
all), we can chose a neighborhood small enough that all intersections are transverse, 
and therefore m will be constant. The only discontinuities occur when a magnetic 
geodesic is tangent to dR, that is exactly at the orbits in Rq. 

We will now look at these discontinuities. If a magnetic geodesic 7 corresponds 
to an orbit b in Ro, it is tangent to R at a point 7(0). It agrees with a magnetic 
geodesic in M, so it never reaches R again and m(b) — 1. By the magnetic convexity 
of M, there is a 6 small enough that each magnetic geodesic in the ball .6,5(7(0)) 
goes through R at most once. Moreover, since the magnetic geodesic is compact, 
we can find an e neighborhood ^(7) such that it only intersects R close to 7(0), 
i.e. JV e (7)ni2 = B 4 (7(0))niJ. 

If the orbit b corresponding to 7 is regular, orbits in a small enough neighbor- 
hood will correspond to nearby magnetic geodesies, completely contained in N e (j). 
These magnetic geodesies will then intersect R only inside £,5(7(0)), and therefore 
at most once. Thus, these orbits will have m equal to or 1. 

We have that B is a 2-dimensional orbifold, and Ro is a continuous curve on 
it. The function m is constant on each connected component of B\Rq, and takes 
values or 1. On Rq, the function m = 1, except maybe at isolated singular orbits. 
Since on regular orbits m = 1, we can say that magnetic geodesies go through R 
at most once, except maybe at a finite number of singular ones. Any singular mag- 
netic geodesies that is tangent must go through R only once. If a singular magnetic 
geodesies cuts dR transversely, we know that m = 1. But the corresponding orbit 
pases through SR exactly m/a times, so a = m = 1 and the geodesic is not singular. 

In the case where SM is non orientable, consider instead its orientable double 
cover SM, and the associated vector field G. Then SM is a compact orientable 
3-dimensional manifold, and G is a smooth vector field that foliates SM by circles. 
We can follow the same arguments with a few modifications. 

The correspondence between orbit on SM and magnetic geodesies is not a 1 — 1 
correspondence, a magnetic geodesic lifts either to an orbit that covers it twice, 
or two disjoint orbits. Nonetheless, will use this correspondence freely, keeping in 
mind this possible duplicity. 

Let SR be the subset of SM that corresponds to the region R, and SdR the 
subset of SM corresponding to vectors tangent to the boundary of R. Let B be 
space of orbits of SM and Ro the set of orbits passing through SdR. The counting 
function m : B — > N can then take value 2, since an orbit that covers a magnetic 
geodesic twice will pas through SR twice. In fact, when m(b) 7^ 0, it will be 1 if the 
magnetic geodesic corresponds to two disjoint orbits, and 2 when it corresponds to 
an orbit that covers it twice. 

Since on regular orbits m = 2 only on orbits that cover a magnetic geodesic twice, 
we can say that magnetic geodesies go through R at most once, except maybe for 
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a finite number of singular ones. Any singular one that is tangent must go through 
R only once, by assumption. If a singular orbit cuts dR transversely, we know that 
m is at most 2. But the orbit pases through SR m/a times, so if it is singular 
a = m = 2 and the geodesic goes through R only once. 

Every magnetic geodesic goes through R at most once, and outside R they agree 
with the magnetic geodesies from M. For the magnetic geodesies to close, they 
have to exit R in the same place and direction, therefore preserving the scattering 
data. □ 

If the region R is simple, we can use this result together with theorem 14. 2 1 to get 
rigidity. For surfaces of constant curvature is easy to see that any circular disk that 
is strictly smaller than one of the orbit circles is a simple domain. Corollary [3] can 
be stated in a more precise way as the following theorem. 

Theorem 5.1. Let M be a surface of constant curvature K , and k > big enough 
that all circles of curvature k are simple. Let rk be the radius of a circle of curvature 
k, and < r < r^ such that r + r^ is smaller than the injectivity radius of M . Let 
R be a compact region contained in the interior of a disk of radius r. Then the 
region R can't be perturbed while keeping all the circles of curvature k closed. 

Consider M with the constant magnetic field that has circles of curvature k as 
magnetic geodesies. If R is contained in a disk D of radius r we can consider any 
perturbation R of R as a perturbation D of D. Since r < the disk D is simple. 
Also, since r + rj, is smaller than the injectivity radius of M, any circle of curvature 
k will go through D at most once. We can then use theorem[2]to show that D and D 
have the same scattering data. Since D is simple, and they have the same scattering 
data, by theorem !4.2l £) is also simple. But in [DPSUl Theorem 7.1] N. Dairbekov, 
P. Paternain, P. Stefanov and G. Uhlmann proved that two 2-dimensional simple 
magnetic systems with the same scattering data are gauge equivalent. 

If M is not compact, consider instead of M a compact quotient that contains 
all the magnetic geodesies that pass through D. This can be achieved since all this 
magnetic geodesies are inside a disk of radius 4r, where r is the radius of a circle 
of curvature k. 
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